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We calculate the diffusion thermopower of the degenerate two-dimensional hole gas in a p— type Si/Sii_a;Ge 3 
lattice mismatched heterostructure at low temperatures and zero magnetic field. The effects of possible 
scatterings, e.g. remote impurity, alloy disorder, interface roughness, deformation potential, and random 
piezoelectric on the hole mobility and the diffusion thermopower are examined. Calculated results are 
well fitted to the experimental data recently reported. In addition, we predict a possibility for the diffusion 
thermopower to change its sign as the SiGe layer thickness changes, the effect has not been discussed yet. 
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1 Introduction 

In a recent experiment, the thermopower of the two-dimensional (2D) hole gas in a p— type Si/Sio.8sGeo.i2 
heterostructure (HS) were reported 1 1 1. At very low temperature, the total thermopower S, approximated 
by the diffusion thermopower S , was explained by a phenomenological model without calculations start- 
ing from the microscopic level JTJ. Since many theoretical studies EJG] of the diffusion thermopower are 
available, it may be interesting if one can use them to examine the reported data. 

Among the scatterings determining diffusion thermopower, impurity doping and interface roughness 
are normally studied G]|3). However, proper considerations for diffusion thermopower due to deformation 
potential and alloy disorder, to the best of our knowledge, are unavailable. Further, some assumptions 
used in the existing studies are actually weak. First, the potential barrier at the HS's interface is always 
assumed to be infinite (2j [3), although it is small, thus changes the transport properties of the HS J3J 
|4). Next, deformation potential scattering limiting the 2D hole mobility is always based on the idea that 
the holes and the electrons undergo the same deformation potential 0. This assumption was indicated 
H to be invalid, and the properly derived deformation scattering was found |4] to strongly limit the 2D 
hole mobility. Consequently, one may interest in studying the diffusion thermopower determined by this 
scattering. Lastly, a study of the diffusion thermopower caused by random piezoelectric field, the recently 
introduced scattering J4], is open. This scattering was found [6| to cause a sign change of the diffusion 
thermopower in an n— type GaAs quantum well as the well thickness varies. Even though piezoelectric 
scattering in SiGe alloy is weak [4J, a study of the diffusion thermopower associated with it may be useful. 

In this paper, we calculate the diffusion thermopower of the 2D hole gas in a p— type Si/Sio.8sGeo.i2 HS 
with finite barrier. The scatterings to be considered are remote impurity, alloy disorder, interface roughness, 
piezoelectric, and the properly derived deformation potential. Calculated diffusion thermopower is given in 
comparison with experiment [ 1 1. We also discuss the possibility for the diffusion thermopower to change 
its sign as the SiGe layer thickness varies. 
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2 Diffusion thermopower 

The total thermopower S is defined by £ — SWT in the presence of a temperature gradient VT and 
an electric field £. At the limit of weak coupling, the total thermopower S of a degenerate 2D gas is 
S = S d + S s , where S d and S s are the diffusion and phonon-drag components, respectively EE] [3]. At 
low temperature T, S d oc T H] El E), while 5 s has a more complex behavior depending on the carrier- 
phonon scattering. For deformation potential scattering, S s oc T 6 [7|, while for piezoelectric scattering, 
S g oc T 4 |8|. In the former case, S s dominates as T < IK, while in the latter, S d dominates below 0.5K. 
Thus, as T < 0.5K, one can consider S d as the total thermopower, which can be measured [ 1 ]. For elastic 
scatterings, the diffusion thermopower S d of a 2D hole gas at zero magnetic field is given by |E1E1 : 



n 2 k 2 B T 
3|e| 



d In a{E) 
d~E 



(1) 

E—Ep 



where |e| is the hole charge, kp Boltzmann constant, Ep = h 2 k F /2m* the Fermi energy, kp = y/2irps 
Fermi wave number, and m* the hole effective mass. In Eq. ([TJ, the a(E) is the conductivity given by 
a(E) = ps{E)e 2 T(E)/m* with t(E) is relaxation time and ps(E) = Em* /nh 2 the 2D hole density 
El El- It is always assumed dEl that t(E) oc E p , so the Eq. Q can be rewritten by Mott formula [2 , 3 1: 

with the scattering parameter p defined by |E1EI : 



P r(E F ) 



dr{E) 
dE 



(3) 

E— Ep 

The hole gas in our HS, as will be seen later, is expected to undergo simultaneously remote impurity, alloy 
disorder, interface roughness, deformation potential, and piezoelectric scatterings. In this case, the total 
relaxation time t(E) is given by the Matthiessens rule |9|: 

1 _ 1 1 1 1 1 

t{E) ~ TRl (E) + tad (E) + t ifr (E) + t dp (E) + t pe (E) ' 

where t R i(E), tad(E), ti FR (E), top(E), and tpe(E) are the relaxation times due to the scatterings 
listed above, respectively. 

It could be seen from the Eqs. ^ and that the total diffusion thermopower S d of a 2D system, 
which undergoes more than one scatterings, is specified not only by the scattering strengths, but also the 
energy dependence of the relaxation time t(E) at E = Ep. In terms of the corresponding autocorrelation 
functions (\U(q)\ 2 ), they are expressed by (4l|6l: 



t(E] (2k)*HEJ » J J 4,2 _ q ,2 t'tff 

where q' — (q' , 9'), E — h 2 q 2 /2m* is the energy corresponding to wave vector q — (q, 6). The dielectric 
function e(q) in Eq. ^ is given at zero temperature within the random phase approximation by [9 ] 

e(q) = 1 + ^F s (q) [1 - G(q)] for q < 2fc F , (6) 

with qTF = 2m*e 2 /e^h 2 is the inverse 2D Thomas-Fermi screening length, ep the dielectric constant of 
the HS, and the function G(q) = q/2^/q 2 + 4/c| in Eq. (j6jl allows for the local field corrections associated 
with the many -body interactions of the 2D hole gas |9|. The form factor Fs(q) is defined by Bl flOl : 

r+ca r+oo 

F s (q)= dz rfz'ICWnC^')! 2 ^- 2 ' 1 , (7) 
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where £(z) and C( z ') are the wave functions representing the two holes having interaction An explicit 
expression for Fs(q) will be given in the Eq. ( 10 1, thus (\U(q)\ 2 ) is all needed for specifying t(E). 



3 Autocorrelation functions for the scattering mechanisms 

Scattering by a random field can be specified by its autocorrelation function in the wave vector space 
{\U(q)\ 2 ) f£|. Here the angular brackets stand for an ensemble average over the fluctuations of the 2D 
Fourier transform of the random scattering field U (q), given by J4]|9l: 

/+oo 
dz\az)\ 2 U(q lZ ). (8) 
-oo 

In our HS, the 2D hole gas is confined by a triangular potential |fl~|[TT) located along the growth direction 
chosen as the z axis with the Si/SiGe interface locates at z — 0. It has been shown lfT2ll that for the finite 
barrier potential Vo, the lowest subband may be very well described by the modified Fang-Howard wave 
function 0H2 [TO): 

, . _ J Ak 1 / 2 6 kz / 2 for z < 0, 

i(Z) ~ { Bk^ikz + c)e- kz / 2 for z > 0, C ' 

in which k and k are half of the wave numbers in the well and barrier, respectively. Further, A, B, c, k, 
and k are the variational parameters to be determined by minimizing the ground state energy connecting to 
Vo lfT2ll . Detailed expressions to be minimized could be found elsewhere l4l l9l [lOl IT2l [T3l . 

Using Eq. (|9|, the form factor Fs(q) in (|7jl can be given in term of t = q/k and a = n/kby iPfl fTOl 

A 4 a „ 9 , 2 + 2c(t + l)+c 2 (t + l) 2 B 4 r . , , , 

Fs{t) = + 2A 2 B 2 a A > J \ L + — 2 c 4 + 4c 3 + 8c 2 + 8c + 4 

w < + a (t + a)(i+l) 3 2ft+l) 3 L v ' 



t (Ac 4 + 12c 3 + 18c 2 + 18c + 9) + i 2 (2c 4 + 4c 3 + 6c 2 + 6c + 3) 



(10) 



Our HS is doped by an impurity sheet of thickness L\ supplying a scattering source JT] [TTJ . This 
impurity sheet follows a spacer layer of thickness Lg grown on the top of the SiGe layer of thickness L 
ifTTI . Thus, the first autocorrelation function we have to specify should be l9l fT2l [131 : 

(WkM?) = — -r- / dzMzi)F 2 (q, si) (11) 

V e Lq ) q + qi Jl+l s 

where F\(q,Zi) is the remote impurity form factor of a 2D impurity sheet located at defined by ll9l fT2l : 

/•OO 

F 1 (q,z i )= / dz\az)\ 2 e-^-^ (12) 



while n\(z) is the impurity density at Zi, which is m as L + L$ < Zi < L + L$ + L\ and zero otherwise. 

Further, it was indicated [ 14] that at high impurity doping level, the impurity distribution is not com- 
pletely random. Due to the Coulomb interaction among the charged impurities during the sample growth, 
the impurity diffusion tends to diminish the probability of large fluctuations of the impurity density [TT3TI - 
The so-called impurity correlation effect can be taken lfT3l[T6l by adding to the autocorrelation function a 
screening-like factor qj (q + qj) in term of the inverse statistical-screening radius qi defined by lfT3l[T6l 

27re 2 ?iiLi 

qi = — tt^- < 13 ) 

In the above expression, Tq ~ lOOOif is the freezing temperature of the impurity system lfl3l[T6l . 
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By plugging Eq. (Nj) into Eq. ( 12 1, and then putting results into Eq. ( 1 1 1, we obtain (\Um(q)\ 2 ) as: 



(\Uw(q) 



f2ne 2 



ni q 



\ e L q J k q + qi [ 2 



1 r 



AB 2 



Qi(§,-) Q 2 (|,5 



i(t,v) 



A 2 a 



y/i {t + a 



2c 



i-t {i-ty {i-tf 



B4 [^(!' s )-^(!' d )]} (14) 

with s = k(L + Lg),d = k(L + Lg + L\) and the auxiliary functions defined by: 



(15) 



Q 2 (M) = 



1 + t 



(c + v) 2 2c(t-3) 2u(i 3 - 3i 2 - i + 3) 4(< 2 - 2t + 3) 

i 2 - 1 + (t 2 - l) 2 + (t 2 - l) 3 " + (t 2 - l) 3 



, (16) 



Qa(t,v) 



e~ 2v t 2 



r 3i 8 - 20t 6 + 66t 4 - 84t 2 + 163 2(c + w)(3i 6 - 17i 4 + 41i 2 - 91) 



(t 2 - l) 6 (t 2 - l) 5 

2(c + «) 2 (3i 4 -14t 2 + 43) 4(c + w) 3 (t 2 -5) 2(c + v)^ 



(< 2 -l) 4 (t 2 -l) 3 (t 2 -l) 2 

The other autocorrelation functions can be taken from Ref. [4|. For alloy disorder, it is given by [4|: 



(17) 



(\U AU (q)\ 2 ) = s(l - aOt&fio 



B 4 6 2 



c%(26) + 4c 3 pi(26) + 6c 2 p 2 (26) + 4cp 3 (2&) + p 4 (2&) 



, (18) 



where a; denotes the Ge content of the SiGe alloy, mai is the alloy potential. The volume occupied by one 
alloy atom is given by ilo = a A1 /8, with <zai the lattice constant of the alloy. The auxiliary functions pi(v) 
(I = — 4) of the dimensionless variables v and b used here are defined by |4| 

a < 1 



J=0 



J! 



(19) 



The autocorrelation function of interface roughness scattering is given by [4| 

PsB 2 



\UlFR(q)\ 2 ) = 



47re 2 x 2 



B 4 



n D c z + 2c + 2 



c 4 + 4c 3 + 8c 2 



Be + 4 



(20) 



where hd is the depletion charge density. (| A 9 | 2 ) is the spectral distribution of the interface profile usually 
assumed 12 ID to be in Gaussian form specified by roughness amplitude A and correlation length A as: 



(|A ? | 2 ) =^A 2 A 2 exp 



q 2 A 2 



(21) 



Deformation potential scattering is a combined effect of lattice mismatch, which gives rise a strain 
field e, and interface roughness ifTTl [T8l . In the previous studies [5|, it has always been assumed that 
the deformation potential experienced by the holes in the valence band is identical to that experienced by 
electrons in the conduction band [ 17 1 with a different coupling constant S. This assumption is, in fact, 
invalid. In particular, while the deformation potential for the electrons in the conduction band is fixed by 
a single component e zz of the strain field, that for the holes in the valence band must be fixed by all three 
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diagonal components of e fl9\ . Thus, the autocorrelation function for the deformation potential scattering 
for the holes needs to be modified to have the following form [4]: 



(\UuAq)\ 2 ) = 



B 2 b 2 ae ]l 
2L 



c 2 p (b + bt) + 2cp! (b + bi) +p 2 (b + bt) 

2 



[b s (K + 1)] 2 (1 + sin 4 9 + cos 4 



/ d s GY 

V4C44/ 



1 



sin 2 29 

~ 4 



(22) 



(|A 9 



in which b = kL, e\\ = (as; — aAi)/aAi 1S the in-plane strain, b s and d s the shear deformation potential 
constants. The anisotropy ratio a = 2c^j(c\\ — c\%) and the elastic constants K — 2cnlc\\, G = 
2(K + l)(cu — C12) are deduced from the elastic stiffness constants of the SiGe alloy en, C12, and C44. 

Similarly, random piezoelectric scattering is a combined effect of lattice mismatch and interface rough- 
ness [4|. The other requirement for this scattering is the piezoelectricity of the strained SiGe layer, which 
has recently been found |20|. Because of the interface roughness, the off-diagonal components of the strain 
field in the SiGe layer become randomly fluctuating HHTT]. Therefore, they induce inside the SiGe layer 
a fluctuating density of bulk like piezoelectric charges supplying a scattering source. The autocorrelation 
function for piezoelectric scattering in our HS has been derived by Ref. [4| as: 



( 37ree 14 Gcce 1 1 sin 29 
V 8e L C44 



(23) 



with ei4 is the piezoelectric constant of the SiGe alloy [20]. The piezoelectric form factor -Fpe(£) 
appearing in the above equation is given by J4] 



A 2 a , 
FpE(t) = — - 1 - e 
t + a v 



-2bt\ 



B 2 b 



2c l t 



At 



3 +c 2 (l - 2bt)p (b + bt) 



Act 

+ (t + i) 2 '• {t + iy 

+2c(l + ct- 2bt)p 1 {b + bt) + (1 + Act - 2bt)p 2 (b + bt) + 2tp a (b + bt) 

-e- 2bt c 2 p (b-bt)+2c Pl {b-bt)+p 2 {b-bt) j. (24) 

The background impurity scattering examined in some previous studies (2J |3] is not considered here. 
The reason is that our HS has only an impurity sheet separated from the well by a spacer of thickness 
L s = 120A, but no intentional background impurity [DQT). Thus, the Eqs. (IT) , |l8} , po} , |22| , and ((23 1 
supply all needed for specifying the relaxation time using the Eq. Q. 



4 Numerical results and discussions 

4. 1 Comparison to experiment 

We now compare calculated S d to some experimental data for the thermopower S (T < 0.5K) of the 
2D hole gas in a p— type Si/Sio.88Geo.i2 HS reported on the Fig. 6 of the Ref. [1 1. This HS, which is 
the sample CVD191 used in the Refs. ifTl fTTIl . is composed by a strained Sio.8sGeo.12 layer of thickness 
L grown on a Si substrate. A spacer of thickness L$ followed by an impurity sheet of thickness L\ is 
placed on the top of the well. At p§ = 2.7 x 10 n cm~ 2 , the hole mobility fj, is 15000cm 2 /Vs, while the 
thermopower S is given [1 1 by the solid squares in the Fig. 1. As T < 0.5K, S is approximately linear, 
reflecting the domination of the diffusion thermopower S d calculated by Eq. ^ using a phenomenological 
expression for p with five fitting parameters but no calculations starting from the microscopic level (TJ. 

Our calculations, on the other hand, are based on experimental parameters using A and A as fitting 
parameters. These parameters are: L = 400A, L s = 120A, L T = 300A, x = 0.12, n T = 3.0 x 10 18 cm~ 3 , 
m* = 0.29m e [DQ3]]. The finite barrier is chosen to be V (x) = 0.74a; = 0.089 eV as in the Ref. ED 
while the alloy disorder potential uai — 0.30 eV as in the Ref. [22 1. Following the Refs. flUEI], the other 
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parameters are: asi — 5.43 A, a Gc — 5.658 A , <zai = 5.455 A, n D = 5.0 x 10 10 cm 2 , cl = 12.192, 
b s = -2.804 eV, d s = -5.240 eV, c n = 16.15 x 10 10 Pa, c 12 = 6.203 x 10 10 Pa, c 44 = 7.821 x 10 10 Pa 
and ei4 = 0.956 x 10 -2 C/m 2 . Detail discussions in choosing parameters are available in Ref. [4|. As a 
result, calculations with A = 1.3 A and A = 97 A give n = 14709 cm 2 /Vs and S d /T is -12.60 ^V/K 2 , 
providing the best fit to the reported data, as seen on the Fig. 1. 

It has been suggested both experimentally 021 US ED and theoretically [13]|24] that A — 100 A while 
A varies from 1 A to 20 A. In fact, the roughness amplitude A = 1.3 A is small comparing to those 
normally used J2][3j|4][5j[6), thus it is necessary to propose an interpretation for this value. It has been 
pointed out [ 17 18 1 that a strained layer, as thinner than a critical thickness [25 1, prefer to relax by buckling 
which gives raise interface roughness. Consequently, A is found to depend strongly on the strained SiGe 
layer thickness [18], the Ge content of the SiGe layer [18 23), and the cap layer thickness 1181 . In our 
HS HHm, the SiGe layer with L = 400 A is thick, the Ge content x — 0.12 is small, and the cap 
layer is thick (Ls + Lj = 420 A). All of them, interestingly, supply a small A lfl8l |23l . Further, there 
exist in the literature many studies supporting A at the same order. Very small roughness amplitudes 
(A ~ lA) of Si/SiGe interfaces have been seen experimentally in the Ref. [23 1. The same A were also 
used theoretically, including A = 2 A in the Ref. 03, and A = 1.78 A in the Ref. Il24l . 



4.2 Scattering mechanisms and the diffusion thermopower 

In order to examine the strengths of the existing scatterings in our HS, the hole mobilities limited by 
separated and combined scatterings are given on the Fig. 2 using the expression fi = eT(Ep)/m* within 
the linear transport theory J9)- The parameters of the Fig. 1 are used here. The Fig. 2 indicates that in 
our HS, piezoelectric scattering is weak and alloy disorder is of minor important on the whole range of 
pg. As ps > 1-5 x 10 u cm~ 2 , deformation potential and interface roughness are the dominant scatterings 
while in the region p§ < 1.5 x 10 11 cm -2 , remote impurity dominates. There are two reasons for the major 
importance of interface roughness, in spite of small A. First, since the barrier is small, the hole wave 
function penetrates deeply into the substrate. Consequently, the hole density at the barrier is finite instead 
of zero for the infinite barrier, thus strengthening interface roughness scattering [4 |. Next, the finite barrier 
suppresses both alloy disorder and deformation potential scatterings, as showed by the Ref. J4). 
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Next, we examine the partial and total S d /T which are plotted vs. p§ on the Fig. 3 with the same 
parameters of Fig. 2. The notations Spp S^ D , Sf FR , Sp)p, Sp E represent the partial diffusion thermopow- 
ers due to remote impurity, alloy disorder, interface roughness, deformation potential, and piezoelectric 
scatterings, respectively. As mentioned above, the total S d is a combination of the partial components, 
weighted by the corresponding scattering strengths. An examination of Fig. 3 reveals that S d changes its 
sign atps — 1-8 x 10 n cm~ 2 , which is smaller than that reported for 2D electron gas J2l[3]. On the whole 
range of p$, Sp E and Sf FR are almost independent of pg. Differently, Sp\ P is very negative and depends 
strongly on p$ as ps < 2.0 x 10 11 cm -2 before becoming independent of p$ as pg > 4.0 x 10 11 cm -2 . 
While S*^ D is very small and may be neglected, S^j is the only one having large positive values. 




1 2 3 4 5 6 7 

1 1 -2 

Hole sheet density p s (10 cm ) 



Fig. 3 Partial and total S ld /T of the Si/Sio.88Geo.i2 HS in Fig. 2 as functions of ps with the parameters taken from 
Fig. 2. RI, AD, IFR, DP, PE are the abbreviations of remote impurity, alloy disorder, interface roughness, deformation 
potential, and piezoelectric scatterings, respectively. 
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Now we turn to another interesting issue: the possibility of changing in sign of S d when the SiGe layer 
thickness L changes. While the change in sign of S d in n— type Si-MOSFETs as the carrier density varies 
has been addressed [2 , 3 1, no discussion on the dependence of S d on L has been given. A recent study 
on the other hand, suggests that the S d of an n— type GaAs quantum well can change its sign as the well 
thickness L changes. This possibility is a consequence of a strong piezoelectric scattering in the quantum 
well made by GaAs material with a large piezoelectric constant ei4. It may be interesting to figure out if 
there is such a possibility in a SiGe alloy with a smaller ex4? 



30 



RI_ 



0.1 * Uoc: 
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1 - t'AD 














Urn 










-20 -EE- 



40 60 80 100 120 
SiGe layer thickness L (A) 



140 



Fig. 4 Plot of partial and total S d /T of a Si/Sio.8sGeo.i2 HS as functions of the SiGe layer thickness L with Ls = 
100A (solid line). Dashed and dotted lines give the total S d /T for L s = 125A and L s = 150A, respectively. The 
hole mobilities limited by separated and combined scatterings of this HS are given on the inset. 



To answer this question, we look for experimentally attainable parameters which allow our HS to exhibit 
a change in sign of S d as L changes. The following parameters are kept: ps = 2.7 x lC^cm -2 , L\ = 
300 A, x = 0.12. The others are chosen for a weaker remote impurity scattering: rij = 1 X 10 cm , 
L s = 100 A, A = 3 A, A = 100 A. The partial and total S d /T are plotted vs. L on the Fig. 4, which 
shows that the total diffusion thermopower changes its sign at the SiGe layer critical thickness Lq ~ 96 A. 

For more information, the inset of the Fig. 4 provides the hole mobilities of the HS limited by separated 
and combined scatterings. It can be seen that while alloy disorder and piezoelectric scatterings are small, 
the others are comparable. Thus, remote impurity, deformation potential and interface roughness are im- 
portant scatterings determining the total diffusion thermopower. Since nothing but S^j 1 (Fig. 3), we 
can adjust Lq by changing the strength of remote impurity scattering. There are several ways to do that, in- 
cluding changing the spacer thickness or impurity density. Indeed, our calculations reveal that Lq depends 
strongly on Ls- For illustrations, the dashed line on the Fig. 4 shows that for Ls — 125 A, Lq ~ 76 A. If 
the spacer is wider, (Ls — 150 A), the critical thickness is even much lower: Lq ~ 55 A (the dotted line). 

5 Conclusion 

In conclusion, we present a theoretical study of the diffusion thermopower S d in a p— type Si/Sii-^Gej; 
lattice mismatched HS at low temperature and zero magnetic field. In the HS, deformation potential, alloy 
disorder, and piezoelectric scatterings are examined in comparing to the conventional scatterings. The 
calculated diffusion thermopower is in good agreement with a recent experiment. Further, S d is found 
to depend strongly on the SiGe layer thickness L, and changes its sign as L across a critical thickness 
Lq. The possible parameters which can affect Lq is also proposed. Deformation potential is a dominant 
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scattering, making an important contribution to S d . On the other hand, piezoelectric is weak while alloy 
disorder has a very small partial S d . Changing the hole density, we find a sign change of S d at a smaller 
hole density comparing to that reported previously for 2D electron gases. 
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